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ABSTRAGT 

Spectrum and chemical composition of cosmic rays accelerated in supernova remnants are 
studied on a basis of kinetic approach. The cosmic ray transport equation with the Bohm 
diffusion coefficient has been numerically solved self-consistently with gas dynamic equations 
for the underlying flow. Gomparison with observational results gives some indications that 
galactic cosmic rays are produced by supernovae shocks which expanded in the relatively low 
temperature partially ionized interstellar medium. 

INTRODUGTION 

Supernova remnants (SNRs) are considered as the main source of cosmic rays (GRs) in the 
Galaxy. Diffusive shock acceleration process is able to convert sufficient amount of the explosion 
energy into GRs and to produce GR spectrum with necessary shape and amplitude (Berezhko et 
ah, 1996). At the same time there is no observational conhrmation that the nuclear component 
of GRs is produced in SNRs. Investigation of GR chemical composition can give additional 
evidence whether SNRs are indeed the source of observed GRs (Berezhko et ah,1995; 1996; 
Ellison et ah, 1997). The most pronounced aspect of the observed GR chemical composition 
is the increasing of relative abundance as function of atomic number A. Here we present the 
results to study under which conditions the theory hts the observed features of GR chemical 
composition. 

METHOD 

The description of SNR evolution and GR acceleration is based on the kinetic approach. It 
includes (Berezhko et ah, 1996) the diffusive transport equation for the GR distribution function 
which is solved self-consistently together with the gas dynamic equations for the underlying flow. 
The model includes an assumption about injection of some fraction 77 <C 1 of superthermal gas 
particles (ions) at the subshock into the acceleration regime. By dehnition gas particles becomes 
GRs when they reach some minimum velocity Vinj during the thermalization. Therefore the 
momentum of injected particles 

PAinj Apinj ( 1 ) 

is proportional to the atomic number A. Here and hereafter variables without subscript A 
corresponds to protons. 

The expected chemical composition of accelerated GRs is sensitive to the dependence of 
the injection rate rjA upon the atomic number. According to the hybrid plasma simulation of 
quasiparallel collisionless shock (Trattner and Scholer, 1993) the relative number of particles 
in the power-law tail of the spectrum is approximately the same for protons and a-particles. 
Therefore we have to take 

9a = 9- (2) 

We assume that supernovae (SN) shock expands into the uniform interstellar medium (ISM) 
with normal chemical composition. Therefore the relative number of accelerated nuclei 
NAinj/Ninj = ( 77 ^ 1 / 77 ) (a^/a) = a^/a is determined by the solar system abundance a a of elements 
with the atomic number A. 


The shape of accelerated CR spectrum essentially depends upon the assumed CR diffusion 
coefficient. We use 

^a{pa) = (3) 

QaG 

which at relativistic energy (pA 3> Amc) coincides with so called Rohm diffusion coefficient. 
Here c is the speed of light, m and e are the mass and charge of protons. The ion charge 
number Qa is assumed to be a function of momentum pA- At the beginning of acceleration, 
that corresponds to small momenta pa ~ PAinjj Qa coincides with the equilibrium ion charge 
state in background ISM Qai- At nonrelativistic energies the acceleration process is extremely 
fast because of the small value of the diffusion coefficient. Therefore we assume the constant 
ion charge number 

Qa{pa) = QAi at Pa < Amc. (4) 

The acceleration process becomes progressively less rapid with increasing CR energy. Therefore 
at the sufficiently high energies pa 3> Amc all kind of ions becomes completely ionized, that 
means 

Qa = Za at pa S> Amc, (5) 

where Za is the nuclear charge number of elements with atomic number A. 

Due to the rigidity-dependent CR leakage from the Galaxy, the expected CR spectrum 
(flux) is connected with the spectrum Ja{pa) produced in SNR (source spectrum) by the relation 

JAobs{pA) oc Ja{pa)t{pa/Za), (6) 

where we assume the power-low rigidity dependent mean CR residence time 

t{R) oc (7) 


RESULTS AND DISCUSSION 

In order to investigate the influence of nonlinear effects produced by the accelerated CR back- 
reaction and ISM parameters on the CR chemical composition we have performed calculations 
for high [f] = 10“^) and low [f] = 10“^) injection rates and for two typical phases of ISM: warm 
(with temperature T ~ lO^K) and hot (T ~ 10®K). The initial ionic charge numbers Qai = 1 
for all elements in warm ISM and Qa slowly increases from QAi = 1 for H and Qa* = 2 for He 
to Qai = 9 for Xe in the hot ISM (Kaplan and Pikelner, 1979). The value of the parameter a 
(see eq.(7)) has been chosen to fit the observed proton’s CR spectrum. 

Calculated hydrogen and helium spectra and CR abundance relative to solar system 
abundance versus atomic mass number (enhancement factor) 


e(A,£fc) 




( 8 ) 


at kinetic energy per nucleon = SAk/A = 3GeV are compared on Figure 1 with experimental 
data. The normalization of calculated spectra was made by htting the observed hydrogen 
spectrum at kinetic energies Sk = 10^ lO^GeV (the relative normalization of hydrogen to 
helium and other nuclei is hxed by the model). In the case of low temperature ISM our theory 
is applicable only for gas-phase ions only for which experimental data are presented on Figure 
lb. We do not consider here the problem of acceleration of refractory elements, which locked 
in grains in the low temperature ISM (this interesting question was considered by Ellison et al, 
1997). 

One can see from Figure 1 that in the case of low injection, when nonlinear modihcation of 
the SN shock by the CR backreaction is negligible the theory in contrast to observation predicts 





Fig. 1: CR spectra in kinetic energy per nucleon (a) and CR enhancement factor versus atomic 
mass number (b) for ISM temperature T = 10® and injection rate rj = 10“® (dash-dotted lines), 
for T = 10"^ and p = 10“® (dashed lines), for T = 10"^ and p = 10“^ (full lines) 


much less efficient production of all kinds of species relative to protons. CR enrichment by heavy 
elements is not essential also for the high injection rate {p = 10“^) in the case of hot ISM. Only 
in the low temperature ISM, with a low degree of ionization for all kind of species, the theory 
predicts essential CR enrichment by heavy elements if the injection rate is relatively high. 

The shape of CR spectrum produced by the modified shock is essentially different at 
relativistic and nonrelativistic energies: J(p) oc p~'^^ at p < me and J(p) oc with 72 < 71 
at p 3> me (Berezhko et ah, 1996). The higher shock modification, the larger difference between 
7 i and 72 . Taking into account the form (3) of the diffusion coefficient, for heavier species we 
have Ja(paM) oc a^(p^/R)“Oi at pa < Qai^c and Ja{pa/^) oc a^(p^/R)“02 at pa 3 > QAi^c 
that gives at relativistic energies the enhancement factor 


e{A) 



( 9 ) 


The ratio A/Za is close to 2 for all elements except hydrogen and a is about 0.7 for three 
considered cases. Therefore for unmodified shock 71 ~ 72 and we have e < 1 that means lower 
heavy element abundance relative to protons. The essential CR enrichment by heavy species 
(e > 1) is expected only under two condition; i) in the case of low temperature ISM with low 
ionization state Qaz — 1 and ii) at relatively high injection rate, which provide the efficient CR 
production and strong SN shock modiheation. It is important to note that according to the 
theory (Ellison et ah, 1997; Trattner and Scholer, 1993; Maikov and Volk, 1995) and experiment 
in the solar wind (Trattner et ah, 1994) the expected injection rate is not less than p = 10“^. 

Only under the above two conditions, as one can see from Figure lb, the theory satis¬ 
factorily agrees with the data. At the same time, as it is seen from Figure la, in the case of 
high injection the theory predicts a much more harder hydrogen spectrum at highest energies 
Ek > lO^^eV than observed. 





Our CR spectra are also essentially harder compared with plane-wave kinetic model pre¬ 
diction (Ellison et ah 1997). The main reason of it is that in our case all phases of SNR 
evolution contribute to the overall CR spectrum (Berezhko et ah 1996). The highest energy 
CRs are produced at the end of free expansion phase, when the SN shock is especially strong 
and produces extremely hard CR spectrum. CRs with intermediate energies Sk = 10^*^ -i- lO^^eV 
are mainly produced during the Sedov phase, when the SN shock is considerably weaker. As 
a result the overall CR spectrum is not of the pure power-law form but has a concave shape. 
In this sense the CR spectra presented by Ellison et al. (1997) correspond to the late Sedov 
phase, therefore they are essentially steeper than ours which include contribution of all phases 
of SNR evolution. 

The possible explanation of the discrepancy between our theory and the experiment can 
be connected either with the acceleration or with the propagation of CRs in the Galaxy. In Erst 
case due to some underestimated physical factors SNR produces a more steeper CR spectrum 
than is predicted by our theory. The second possibility is that the actual escape of CRs from 
the Galaxy is more complicated than described by simple power-law Eq.(7). 
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